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We make several remarks on the B-JIMWLK hierarchy. First, we present a simple and instructive 
derivation of this equation by considering an arbitrary projectile wave function with small number 
of valence gluons. We also generalize the equation by including corrections which incorporate effects 
of high density in the projectile wave function. Second, we systematically derive the dipole model 
approximation to the hierarchy. We show that in the dipole approximation the hierarchy has a 
simplifying property that allows its solution by solving the classical equation followed by averaging 
over the ensemble of initial conditions. 
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Recently there has been a lot of interest in the Balitsky- JIMWLK high energy evolution 0, 0i Q . On one hand 
^vq [ a lot of analytical and numerical work has been done to understand the mean field version of the hierarchy due to 
Kovchegov Q - see for example 0, @ . On the other hand recently first steps has been taken to extend this hierarchy 
in order to include corrections away from the large density limit which are thought to be important for generating 
Pomeron loops 

Still some analytical aspects of the hierarchy itself are not well understood. For example, the original derivation 
due to Balitsky is for a probe containing a fixed small number of partons. A general derivation for a small probe with 
■ an arbitrary wave function has not been given. Such a derivation would be useful not only to polish the surface, but 
also to help understand how to go beyond the small probe approximation. Another aspect that has been discussed in 
the literature 01, IToj but not in a completely systematic way is how to get from the general JIMWLK equation to the 
dipole model [ill], or large N c limit. In particular it is interesting to see whether the dipole limit allows one to probe 
numerically some of the fluctuations which are not included in the mean field approximation of Q , in a reasonably 
Q H ' accessible way. 

This paper contains two sketches which deal with these issues. In Section 2 we present the derivation of the JIMWLK 
qj- equation for a general projectile wave function. We go beyond the standard JIMWLK expression, by including terms 
which arise when the gluon density in the projectile is not small. Our derivation is purely algebraic and very simple, 
r* ■ We hope this approach will be useful for casting the whole problem in a simpler framework and thus facilitate future 
. £h ! progress. In Section 3 we discuss the dipole model limit of the evolution. We derive explicitly the dipole form of the 
evolution and also the 1/N C corrections to it, starting directly from the JIMWLK equation. We make an observation 
that the structure of the dipole evolution is such that it allows the solution of the whole system by solving Kovchegov's 
equation for a single dipole <S-matrix with subsequent averaging over the ensemble of initial conditions. We suggest 
some initial distributions which directly probe fluctuations in the target wave function. 



II. THE S - MATRIX OF A COMPLEX PROJECTILE. 



In this section we derive a compact and intuitive expression for the evolution of the S-matrix of a fast projectile 
which scatters on the fields of some hadronic target. The equation we obtain reduces to the B-JIMWLK 2 equation 
when the projectile is small. For arbitrary projectile we also keep certain corrections due to finite density of the 
projectile partons. 

Although the following considerations can be applied straightforwardly to any colored particles, for simplicity in 
the following we consider projectiles which contain only gluons. We will use the formalism of [T^ . 
Consider an arbitrary projectile whose wave function in the gluon Fock space can be written as 

|*)=*K ta (x)]|0) (2.1) 

The gluon creation operator depends on the transverse coordinate, and also on the longitudinal momentum k + . The 
projectile is assumed to have large energy. The gluon operators which enter ea. (|2.1|> therefore all have longitudinal 
momenta above some cutoff A. We will refer to these degrees of freedom as "valence". Henceforth we omit the 
dependence on longitudinal momentum in our expressions, as the momentum enters only as a spectator variable and 
only determines the total phase space available for the evolution. The wave function in ea. (|2.1fl is arbitrary. In 
particular it does not have to have a fixed number of gluon, but can be a superposition of states with different gluon 
numbers. 
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The projectile scatters on a target, which is described by some distribution of target chromoelectric fields. We 
denote an S-matrix of a single gluon by 



S ab (x) = (0\a?(x)Sa\ b (x)\0) 



(2.2) 



where S is the second quantized S-matrix operator of the field theory. The one gluon S-matrix S a b does not depend 
on the polarization of the gluon and is diagonal in the transverse coordinate x. It does depend on the color fields of 
the target. One can think about S a b{x) as the eikonal S'-matrix, although our derivation in no way depends on this. 
Our only assumption is that the gluons in the projectile scatter independently of each other. Given eas. H2.H2.2(l . the 
formal expression for the S - matrix of the projectile can be written as 



E[5] = (* 



ma?(x)]*[S ab (x)at b (x 



(2.3) 



As discussed at length in|8(, the matrices S(x) are quantum operators on the Hilbert space of the target, and thus 
do not commute with each. To take properly into account this noncommutativity we were led in [8j to consider S as 
a function of additional coordinate x~ and to include the Wess-Zumino term in the weight function for calculating 
averages of functions of S. However for large targets this is not necessary, as the commutator of S is suppressed 
at large density. In the present paper for simplicity we assume that the target is large, and therefore treat S(x) as 
classical c-number fields. The distribution of the target fields is defined by a probability density VFfS], so that 



J dSW[S] = 1, (E) = J dS E[5] W[S] 



(2.4) 



Thus the forward S'-matrix element for the projectile |^) on the target W is given by S = (E). 

To calculate the evolution of the S-matrix we write down the evolution of the projectile wave function to first order 
in the change of rapidity SY. This is given by ^| 



\V(Y + SY)) = 



j-SY I d 2 z(bt(z,[p})bt(z),[p}) 



+ i / d 2 zb1{zM) 



dk+ 



4 a (k + ,z)j\^(Y)) 
(2.5) 

Here the field b[p] depends only on the valence degrees of freedom. It is determined as the solution of the "classical" 
equation of motion 



(l-SY)A V^|fc + I 1/2 



The valence charge density operator is defined as 

P a {z) 



d l b^z)+gr bc b b l {zmz)^gp a (z) 
e ij [d i W-d j bt+gf abc b b M = 



dk + af{k + ,z)T b a c a' i (k + 1 z) 



(2.6) 



(2.7) 



where T h a c = if abc is the SU(N) generator in the adjoint representation. 

The expression ea. (12.51) is not complete when the charge density is large 01 . However it does contain all leading 
terms at small p and in addition certain higher order corrections, which arise due to the fact that b solves the equation 
eq. (12. nonperturbatively. 

The wave function of the projectile after scattering on the target field is 



2tt 



-SY 



d 2 z (bUz t [Sp[) bf(z,[Sp})) 

A 

i I d 2 zb«(z,[Sp}) 



dk' 



T7I ^"(z) a f(fc+,z) *[S ab (x)at b (x)]\0) 



l(l-8Y)A V^\ k+ \ 1/[ 

Here the field b[Sp] solves the equation ea. l|2.6fl but with the rotated source term 

dib1{z,[S P }) + 9 r bc b b (z,[Sp}) b?(z,[Sp]) = gS ab (z) P b (z) 
e« [dib^Sp] - d 3 b«{Sp] + gf abc b b [Sp}b c J [Sp}] - 



(2.8) 



(2.9) 
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Since the initial projectile wave function does not contain soft gluons which appear in ea. 12.5fl . the calculation of 
the averages of the soft gluon operators is straightforward. The S-matrix of the boosted projectile is therefore (to 
first order in SY): 

Zy+6y[S] = Z Y [S] - (2.10) 
' )Y (Q| **[«?] / d 2 z [b?(z,[p])bUz,[p]) + b?(z,[Sp])bUz,[Sp]) - 2b«(z,[p])S ba (z)b b (z,[S P })} nS ab al b }\0) 



2vr 

We now want to represent this expression in terms of functional derivatives of S. To this end we note that 

p a (xMSal}\0) = T b a c 4\x)-^nSal}\Q) = -tr\s(x)T a -^—)nSal}\0) ; 

Sa] c (x) I SSt{x)) 



p b (y)p a (xMSal}\0) = \ T^ d at c (y)-^—\ \t^{x)-^— \*[S4]\0) 

Sal (y) I I Sal J (x) 



= *\ 8 ^mi)r\ 8 ^mr)r l8a 

[S(x)p(x)] a nSa\}\0) = -tr|T°5(!c)^y}*[SaJ]|0) ; 
p b (y)[S(x)p(x)} a nSal}\0) = tr Sr a S(x)-^^\ tr j,%)T b ^L^j *[5o+]|0) etc. (2.11) 

The structure of these relations is quite interesting. The action of the charge density p on the wave function is 
equivalent to the action of the generator of the right rotation SUr(N) on the single gluon S'-matrix, while the action 
of Sp is equivalent to the action of the left rotation SUl{N). In addition the ordering of the action of the SU(N) 
generators on S is opposite to the ordering of p acting on a. 

Using eas. (|2. 1 lfl we can rewrite ea. H2.lO0 as 

^E[S] = X £[S] (2.12) 



bf(z, [J R ]) b?(z, [J R ]) + b?(z, [J L ]) 6?(z, [J L ]) - 2 S ba (z) bf(z, [Jr]) bl(z, [J L ]) (2.13) 



Here 



r R ( x ) = -trl [ S(x)T a J ^Y ja ix) = _ tx ^ T a S(x) _A_Y ja [x) = [S ( x)Mx)] a } (2 . 14) 

and the tildas over b indicate that the ordering of J's in b[J] is opposite to the ordering of p in the solution of ea. (|2.6|) 
- b[p] . Using ea. (|2.4|l this can be recast as the evolution equation for the target probability density W. To do this 
first of all we write ea. (|2.4(l as 

<S)y = j dS ^Y-Y [S] W Yo [S] (2.15) 

where we think of the target as being evolved to rapidity Yq, while the projectile takes on itself the rest of the rapidity 
of the process [1^ . The derivative of £ with respect to Y can be replaced by the derivative with respect to Yq , and 
the physical requirement of Lorentz invariance of the scattering amplitude allows us to differentiate W rather than 
S. Subsequent integration by parts of the functional derivatives with respect to S gives the evolution equation for W 

W[S] = x W[S] (2.16) 



1 

27 



d z z \b?(z, [J R ]) b?(z, [J R ]) + bUz, [Jl]) b?(z, [J L ]) - 2 b b {z, [JlM(z, [Jr]) S ba (z)] (2.17) 



where now the original ordering is restored in b, due to the integration by parts. 
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Let us take the limit of our formulae for small projectile. This amounts to approximating b by the first order term 
in g 

™ = (2 - 18) 
The evolution operator \ reduces to the kernel ^- JIMWLK 

JIMWLK = _^sf (* " X)i(z - y)i r ( s{x)T a_$\ f g ( )T » *\ (2 . 19) 

In this approximation the evolution equation 1(2.16(1 for the probability distribution W becomes the JIMWLK equation 

A W[S] = X JIMWLK W[S] . (2.20) 

To see that this is indeed the JIMWLK equation we recall that a simple form of the latter is written in terms of 
functional derivatives with respect to a(a;,y)[l3 defined as 



S(x) = Pexp^iJ dx-a a {x.x-)T a ^ (2.21) 



with P denoting the path ordering. When acting on any function of S, the derivative with respect to a(x, Y) is 
identical to the action of Jr (x) : 

6Sab (*) 5 men ,rc^™i«6 5 



6*W) m = S^)WHx) F[S] - ^ T ^5S^) F ^ - - U ^ )F[S] (2 - 22) 
Thus we can write ea . 112.20(1 in the familiar form 

w w & - &J m llZZrt-'Jrtoj'.Y) [1 + s,(a;,S(! ' , - sHx)m - sU ' )Slvr W) w[s] (2 - 23) 

We stress that the introduction of the additional coordinate x~ and the path ordering in cq. 12.21fl serves only 
to mimic the action iJr on a function of S as in ea. l(2.22|l . In the framework of B- JIMWLK 2 equation it has no 
independent meaning, since none of the physical observables depend on x~ . 

This completes our derivation of the B- JIMLWK 2 equation. We note again, that eq. I|2.17|l is more general than the 
original B- JIMWLK 2 . Expansion of the solution for the field b in powers of p furnishes corrections to this equation 
due to the high density of gluons in the projectile. 

X - X JIMWLK + 0((a sP 2 ) 2 ) 
We now turn to the discussion of the dipole model limit. 



III. THE DIPOLE MODEL LIMIT. 

In the large N c limit the high energy evolution can be recast in the form of the evolution of Mueller's dipoles[Tl||. 
In this limit one thinks of the gluon as built out of a fundamental and antifundamental "constituents" q and q. The 
wave function then is assumed to have the dipole structure where the color charge is neutralized on pairs of points in 
the transverse plane, namely ^[q a (x)q a (y)]. The wave function of such a superposition of dipoles after the scattering 
event is ^[q a (x)U af3 (x,y)q l3 (y)), where U af3 (x,y) = {Sp(x)Sp(y)} al3 and Sp is the fundamental representation of 
the single gluon scattering matrix. In the large N c limit in calculating £ we neglect all the terms with contractions 
between g's not belonging to the same pair, and therefore 



£ = £[s] ; W = W[s] , 



(3.1) 
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where s(x,y) — jj-Tr[U(x,y)}. Given the form of the B-JIMWLK 2 equation we can derive the evolution of the 
function E, or alternatively W and we can also see how the dipole evolution breaks down due to 1/N C corrections. In 
line with the dipole model ideology we can write down the action of the left and right rotation generators as 

Ju(y)m = ±r tv[S\u)S(v)r a ] [S(v -y)- 5{u - y)\ (3.2) 
N c os{u,v) 



JKvMs] = i tr[S\u)r a S(v)} [S(v - y) - S(u - y)] 

l\ c os(u,v) 

J a R (x)J a R (y)j:[s} = ±r tr[SHu)S(v)r a T a } [6(v - y) - S(u - y)] [S(v - x) - S(u - x)] 



iV c ' Ss(u, v) 

<5 2 E 



+ -L tr[S\u)S(v)r a } tr[St(p)S(r)r°] [S(v - y) - S(u - y)\ [5(p - x) - S(r - x)} 

iVjf os(u,v)os{p,r) 

r L (x)J a L (y)m = i- tr[S^(u)S(v)r a T a ] [S(v -y)- S(u - y)] [S(v - x) - S(u - x)] ^ 



N c " Ss(u,v) 

<5 2 E 



+ -L tT[SHu)r a S(v)} tv[SHp)T a S(r)} [5(v - y) - S(u ~ y)} [S(p - x) - S(r - x)] 

iVjf os(u,v)ds{p,r) 

J b L (x)J a R (y)ns] = 1. tr[SHu)T b S(v)r a } [6(v - y) - S(u - y)] [6(v - x) - S(u x)\ ^ 



N c 1 ' 8s{u, v) 

+ -L tiiSHu)S(v)r a ] tr[5t (p)r b S(r)} [6(v - y) - S(u - y)] [8(p - x) - S(r - x)} 



iV 2 1 ' ' ' 5s(u,v)Ss(p,r) 

Here r a are the generators of SU (N) group in the fundamental representation. We have also dropped the subscript 
F for brevity, but all the matrices in ea. ()^.2j) are in the fundamental representation. To simplify this further we need 
to use the completeness relation 



a a 1 



1 * c 



(3.3) 



and the identity 

Sf(z) = 2Tr [r a S F (z)r b S F (z^] (3.4) 
Working out the color algebra we find that the evolution operator -^ JIMWLK can b e decomposed into two pieces 

x jimwlk [s] = Xdipole[s] + -L Xcc [ s , S ] (3.5) 

The dipole evolution operator Xdipoie involves dipole degrees of freedom only and it is proportional to a single derivative 
with respect to s 

Xdi P oie[s] = 7T f ( } X \if > _ ^ 2 [- S ( X >V) + s(x,z)s(y,z) ] 6 (3.6) 

* n Jx,y,z \ x z ) \ z V) os \ x iV) 

The dipole evolution operator in the form l|3.6[l was first obtained in Ref. |16| from a probabilistic approach to the 
dipole evolution (see also Refs.0,0]). When acting on projectile's E, the operators s and ^ could be viewed as the 
dipole creation and annihilation operators respectively 16] . On the other hand, when acting on W it is more natural 
to write Xdipoie as 

Xdipole[s] = ^ / " %rr^ \2 X f \ [~ S ( X > V) + S ( X > Z ) S (V' Z )] ( 3 - 7 ) 

2tt Jx,v,z {x-zy(z ~ yy ds{x,y) 

and view s as the target dipole annihilation operator while ^ as a target creation operator. 
The operator which accounts for color correlations has the form 

X«= oESr I K(u,v,p,r,z) [Tr(SHu)S(v)SHr)S(p)) - Tr(S^ (u)S(v)S\z) S(r)S*(p)S(z)j\ f 

(3.8) 
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and 

(z - v)i (z - u)i 



K(u,v,r,p,z) 



(z — v) 2 (z — u) 2 



(z-p)i (z-r)i 



(z — p) 2 (z — r)' 2 



(3.9) 



As is obvious from eg. l|3.8fl . the 1/N C corrections cannot be formally represented in terms of dipole degrees of freedom 
only and require introduction of higher order color correlations (see also Ref. |l7)L Thus strictly speaking even if 
at initial rapidity W depends only on s, at rapidities of order N c this will not be the case any more. It is however 
unclear how important is this effect in practical terms. 

Restricting ourselves to the large N c limit we obtain the following evolution equations for the S-matrix E and the 
target distribution function W 

^ = XdipoleMns], dW ^ y [S]) = XdlP ole[s]W(Y,[ S }). (3.10) 

Ea. H3.1Ufl exhibits a curious quantum mechanical-like structure. The evolution in rapidity is analogous to time 
evolution, while the evolution operator %\ is analogous to the Hamiltonian. We can then think of ea. H3.10fl as a 
Schrodinger equation for the wave function W . The peculiarity of this quantum mechanics is that the evolution 
operator ixdipoie is explicitly not Hermitian. This is of course because in our case the probability density is given by 
the "wave function" W itself, and not by its square. The nonhermiticity of Xdipoie in the usual quantum mechanical 
sense is mandatory to preserve the total probability J DsW[s] = 1. Nevertheless, also in the present case one can 
sensibly talk about the Hilbert space of the real functionals of s, and define operators on this Hilbert space. The time 
evolution of these operators is given by the usual Heisenberg equation^] 



dO 
dY 



x,o 



(3.11) 



All that has been said so far applies equally to the dipole model and to the generalized B- JIMWLK 2 equation discussed 
in the previous section. The peculiarity of the dipole system is in the fact that the "Hamiltonian" Xdipoie is linear 
in the momentum conjugate to s. Thus the Heisenberg equation for the operator s is purely classical - the time 
derivative of s commutes with s itself: 

ds(x.y) f , 9 (x — y) 2 . „, . /„,„n 

- ' [s(x,y) - s(x,z)s(y,z)} (3.12) 



dY "J (x-z) 2 (y-z) 2 

Additionally, the observable we are interested is the S-matrix, which depends only on s and not on its conjugate 
momentum. Thus the functional equation for £ can be replaced by an ordinary differential equation for s. The S 
matrix at arbitrary time Y has the form 

S(Y) = [ Ds E[s(Y)] W[s\. (3.13) 



Where s{Y) is the solution of the classical equation ea. (|3.12fi which satisfies the initial condition s(Y = 0) = s. 

Thus we see that there exists an economical way of solving the whole dipole hierarchy without relying on the 
mean field approximation of Kovchegov. One solves the classical equation for the dipole S'-matrix s, calculates £ on 
this solution at arbitrary rapidity Y, and subsequently averages over the ensemble of initial conditions W[s]. This 
numerical program is feasible and is currently being pursued. Compared to the numerical solution of the full JIMWLK 
equation (Ref. 19j) this procedure is much simpler and easier for realization as it is confined to the large N c limit. 
Our goals would be to disentangle the role of the target correlations, which can be probed within the large N c limit 
from the effects due to the 1/N C corrections. The latter can be reliably estimated only by comparing solutions of 
eci. H3.13f) with solutions of the JIMWLK equation, provided identical initial conditions are used. While our numerical 
results will be reported elsewhere, in the rest of this paper we discuss some reasonable choices for the initial ensemble 
W[s}. 

The choice of the initial distribution W[s] is of course tantamount to specifying the complete wave function of the 
target (or at least its large N c limit). As such it is a nonperturbative problem which can not be solved by our present 
means. Nevertheless, as a result of the evolution we expect the S matrix to have fairly robust and universal structure. 
It therefore makes sense to try to model some broad features of the initial distribution which will be reflected in 
the preasymptotic behavior of S. One thing we can probe with the distribution rather than with a single initial 
condition, is the importance of the fluctuations of the target fields for the scattering matrix at large Y. There are two 
basic types of fluctuations one can consider. First, the fields in the hadron are correlated only over finite transverse 
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distance. Thus the scattering amplitudes of dipolcs which scatter at the distance larger than this correlation length 
from each other, must be uncorrelated. To model these fluctuations the initial ensemble must contain configurations 
which are impact parameter dependent. The other type of fluctuations is due to the fact that even at the same point 
in the transverse plane, the fields are not constant in time. The gluon density for example is bound to fluctuate as 
any other quantum mechanical variable. These fluctuations can be modeled within ensemble which is translationally 
invariant configuration by configuration. The latter is much easier to handle numerically. As a simple model of this 
kind consider the following form for the dipole scattering amplitude 

s(x,y) = e - (x - y)2G (3.14) 

This is the McLerran-Venugopalan form of the dipole scattering amplitude pol | , where G has the meaning of the gluon 
density. We can allow for the fluctuations in G by simply integrating over G with the weight of the form: 

W[G] = Woe' (G ~A n) (3.15) 

This distribution allows the gluon density to fluctuate around some central value Go with the standard deviation A. 
The proper normalization of the distribution with this ansatz is 

f r°° (g — g i 2 

1 = J DsW[s] = Wo J dGe a- 2 - (3.16) 

which determines the coefficient Wo as 

Wo = 2/(V5FA(l + Erf[G /V&\) (3.17) 
This simple initial distribution is summarized in the following expression for the S'-matrix 

S(Y) = W J DsdG8{s{x,y) - e -(*-») 2G ) e - (G ~A Go)2 E[s(Y)}. (3.18) 

with s(Y) given by the solution of the classical equation ea. (|3.12(l . For simplicity one could consider a projectile 
consisting of a single dipole 

E[a] = s(x,y) (3.19) 
or for example the perturbative wave function of a virtual photon 

- / d 2 (x - y) Pi* [(x - yf, Q 2 } s(x, y) (3.20) 



Here P 1 stands for the well known probability to find a dipole of the size \x — y\ in the wavefunction of the photon 
with virtuality Q. 

By varying the width of the distribution A one can see how important are the local fluctuations of the gluon 
density. Obvious modifications of the (x — y) dependence, like the perturbative G — > Gln[(ir — y) 2 A 2 ] or anomalous 
dimension G — > G[(x — y) 2 A 2 ] 7 are also straightforward to implement. This distribution is translationally symmetric 
configuration by configuration, and thus does not take into account the finiteness of the correlation length A of the 
target fields. One can however think of it as a valid model for small size dipoles, such that \x — y\ < A. 

A more realistic distribution that depends on the impact parameter can also be constructed. To do this we must 
consider both G and Go as b = [x + y)/2 dependent. For heavy nuclei, the impact parameter dependence of Go can be 
considered as given by the Wood-Saxon parameterization. We should also allow the variance of G to be uncorrelated 
at distances larger than some scale A. In other words we would like to have 

(A(ft)A(O)) = A 2 e- fe / A (3.21) 

this can be achieved simply by distributing A with a Gaussain measure 

expi-i / A(b)D(b-b')A(b')\ (3.22) 
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b.b' 



Where D(x) is the operator inverse to the correlator ea. l|3.21|l . The precise form of D is not particularly important, 
as long as its inverse decreases beyond some scale A and approaches a constant value Aq at zero separation. With the 
appropriate normalization constant 

N(b) = J DAW^ib) ejtpj-ijT [A(6)£>(6-6')A(6')]} (3-23) 



the expression for the S'-matrix becomes 

S(Y) = J DsoDGDA JV -1 [6] x (3.24) 

/ , ,2^,,,\ (G(i,)-G (t,)) 2 j-i f [A(6)£>(6-6')A(6')]| 

S (s (x,y) - e-( x -rt G ^J e STO e l 2ji >."' 1 J T,[s(Y)} . 

The numerical implementation of this ^-dependent measure would be very interesting, but is much less straightforward. 
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